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We demonstrate that a tight transverse trap with the local frequency, ui±, gradually varying in 
the longitudinal direction (x) induces an effective potential for one-dimensional solitons in a self- 
attractive Bose-Einstein condensate. An analytical approximation for this potential is derived by 
means of a variational method. In the lowest approximation, the potential is N(S + l)ui±(x), with N 
the soliton's norm (number of atoms), and S its intrinsic vorticity (if any). The results can be used 
to devise nonuniform traps helping to control the longitudinal dynamics of the solitons. Numerical 
verification of the analytical predictions will be presented elsewhere. 

To be published in Physics Letters A 
I. INTRODUCTION 

A Bose-Einstein condensate (BEC) with attractive interactions between atoms (negative scattering length) may be 
stable in an external trap if the number of atoms in the condensate is below a collapse threshold [l| . In this case, the 
BEC in a nearly one-dimensional (ID) "cigar-shaped" trap, which features tight confinement in the transverse plane 
and a weak potential along the longitudinal axis, can form stable matter- wave packets in the form of bright solitons. 
A single soliton |2| and multi-soliton complexes |3j were created in the 7 Li condensate loaded into a strongly elongated 
optical trap. More recently, solitons whose shape is nearly three-dimensional (3D), were observed in a post-collapse 
state in a condensate of 85 Rb atoms (as well as in 7 Li, in this experiment the sign of the interatomic interactions 
was switched into attractive by means of the Feshbach resonance). 

A fundamental equation which provides for a very accurate description of the dynamics of a rarefied quantum gas 
of Boson atoms in the mean-field approximation is the 3D Gross-Pitaevskii (GP) equation [1J. In a normalized form, 
the equation is 
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where ip is the single-atom wave function, x and y, z are the longitudinal (axial) and transverse coordinates (the 
transverse Laplacian acts on y and z), u>± is the frequency accounting for the tight transverse confinement, U(x) 
a loose axial potential, and g a scaled nonlinearity constant (in the case of self-attraction, g is negative). 

In the experiment, the transverse trapping potential may be axially nonuniform, which corresponds to lu±_ = lu±_ (x) in 
Eq. It may also depend on time, hence lu±_ = lj±(x, t), in the most general case. Actually, the axial nonuniformity 
is an unavoidable feature of any experimental setup, and, on the other hand, specially designed nonuniformity may 
be used as an additional tool for the control of dynamics of trapped solitons. The objective of this Letter is to derive 
an effective longitudinal potential induced (in addition to the explicitly present potential, U(x) in Eq. JQ) by the 
x-dependence of cu±, and the corresponding equation of motion for axial solitons. 

If the transverse potential is much stronger than the longitudinal one, it is natural to reduce the full 3D GP equation, 
0, to an effective ID equation. Different approaches were proposed to achieve this purpose under different 

conditions. In most cases, the reduction is based on assuming a factorized ansatz for the wave function 

ip(x,y,z,t) = exp (- I \ r X ] %\ > ( 2 ) 
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where r 2 = y 2 + z 2 , the transverse width a and local amplitude / being slowly varying functions of x and t. Using 
the Lagrangian representation of Eq. pjl. effective ID equations for a and / can be derived as variational equations. 
First, a is eliminated in favor of /, 



a 2 = to- 1 



yjl + |/| 2 , (3) 
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and then a closed-form ID GP equation with nonpolynomial nonlinearity is derived for / 
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Equation Q with g < admits stable solitary- wave solutions, which were analyzed in Ref. |(j too. 

In the case of the weak nonlinearity, (g/2irui±) |/| 2 <C 1, the nonlinear term in Eq. (£Q) may be expanded in 
powers of |/| 2 , which leads to the cubic GP equation with an additional quintic term that corresponds to higher-order 
self-attraction, 
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In a more direct form, a ID equation with the cubic-quintic (CQ) nonlinearity of this type was derived and employed 
in Refs. A formally similar equation of the CQ type (but with a self-defocusing quintic term) was considered 

earlier in works aiming to take into regard three-body collisions in BEC |l0| ) . Despite the possibility of collapse 
induced by the self- focusing quintic term in the ID setting, the CQ equation has a family of exact soliton solutions, 
which are stable against small perturbations (these new solutions were obtained as an analytical continuation of 
well-known soliton solutions [lj] to the CQ equation with the self-defocusing quintic term). 

The approach based on the factorized ansatz @ was generalized in Ref. [8j to describe configurations with intrinsic 
vorticity, with the ansatz replaced by 

1>(x, y, z, t) = r s exp (- - + iSO^j f s (x,t), (6) 

where 9 is the angular coordinate in the (y, z) plane, and S = 1, 2, ... is the integer vorticity (which naturally carries 
with itself the pre-exponential factor r s ). Actually, S plays the role of the "spin" of effectively ID solitons generated 
by the ansatz. This way, an equation for fs{x,t) similar to Eq. can be derived. 

The Letter is organized as follows. In the next section, using the variational approximation |13| . we develop a 
framework for the analysis of the soliton's dynamics in the model with the axially nonuniform trapping, lu± = to±(x) 
(the variational approach employs the Hamiltonian, rather than Lagrangian). In Section 3, we derive a final result, 
viz. , an effective axial potential for the soliton induced by the x-dependence of u± , and the corresponding equation of 
motion for the soliton. The analytical results are obtained under the natural assumption that a scale of the variation 
of (jJ±(x) is much larger than the size of the soliton. The paper is concluded by Section 4. 

II. VARIATIONAL ANALYSIS 

Aiming to derive an effective potential and equation of motion for solitons in the case of lo± = u)j_(x), we skip 
the derivation of an effective ID equation, and instead adopt a 3D ansatz for the soliton which, in the general case, 
includes the intrinsic vorticity (cf. Eq. ©): 

^ so i = Ar s exp (--^ + iS6^j sech e** (7) 

with amplitude A, longitudinal width W, central coordinate £ and phase 4>, in addition to the transverse width, a, 
and spin, S = 0, 1, 2, that were defined above. The norm of this ansatz is 
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To derive an effective potential for the soliton as a function of the coordinate £, we use the Hamiltonian of three- 
dimensional equation 
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The substitution of ansatz (JJJ in H and straightforward integrations yield 
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This result was obtained under the condition that, as said above, tu±(x) varies on a scale which is much longer than 
the soliton's width, W, the same being assumed about the axial potential, U(x). For this reason, uj\(x) and U(x) in 
Hamiltonian l|10[l are taken at x = £. 

The soliton's peak density, A 2 , can be eliminated, in Eq. i|l(J|) . in favor of the conserved norm, N, using Eq. JSJ. 
This leads to 



N 
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_ (25)! 
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In this expression, the transverse and longitudinal widths of the soliton, er and W, are free parameters. The soliton 
chooses their values by minimizing the Hamiltonian, which leads to conditions 

dH dH 

dW~d(a*)~°- {3) 
The first condition takes a simple form, making it possible to eliminate the longitudinal width, 

W=~~, (14) 
C s gN 

which, obviously, makes sense only for g < 0. Indeed, bright solitons may only exist in the BEC with the attractive 
nonlinearity, corresponding to negative g (unless a periodic optical-lattice potential is present, that can support gap 
solitons at g > [lj,[l5|]). The substitution of W from expression ltTil) in Eq. ifTTl) yields 



N 
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Further, the second condition from Eq. i|13[l amounts to a cubic equation for o~ 2 : 
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It is easy to see that Eq. ()16|l has physical (positive) solutions for a 2 if the normalized number of atoms is smaller 
than a critical value, 

V3C^g J uj± 

The existence of largest iV beyond which stationary solutions do not exist reflects an obvious fact that the underlying 
axisymmetric 3D GP equation gives rise to collapse if N is too large 0] . Equations l(T7|) and (|T2)l predict increase 
of 7V cr (for given lu±) in the lowest vortex state (S = 1), in comparison with its zero-vorticity counterpart, by a factor 
of 2^2, which is a known effect too, see Ref. ^3 an d references therein. 

Alternatively, Eq. (|17l) demonstrates that, for given N, physical solutions exist if the transverse confinement is not 
too strong (otherwise, over-squeezing of the condensate will lead to the collapse): 

- <( -»-" s ^If' (18) 

If there is a region where lj±(x) exceeds (^_L) max , a moving soliton crossing into this region will blow up due to the 
collapse. 
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Further analysis of Eq. (|16fl demonstrates that there are two solutions for a 2 , one decreasing with u± , and the other 
one increasing. The former behavior is natural (stronger squeeze leads to a smaller transverse size of the condensate), 
while the latter one is not. To all appearance, the latter solution is unstable. Therefore, we only take into regard the 
solutions with a(uj±) a decreasing function. In particular, it is easy to demonstrate that this solution yields, for given 
N, the transverse width in the range of 

C s \g\N _ 

— , = cr min < a < oo . 

ttVI + 5 

Here, a mlu corresponds to lo± — (wi) m31[ in Eq. (|18|l . 



III. THE EFFECTIVE POTENTIAL AND EQUATION OF MOTION FOR THE SOLITON 



If a 2 can be found from Eq. then W must be taken as per Eq. (|14|l . and both substituted in Hamiltonian 

. An explicit result can be obtained in the case of a relatively weak nonlinearity, 

uo±(gN) 2 < 
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when a relevant solution to Eq. (|16|l . calculated to first two orders of the perturbative expansion, is 

- 1 {CsgNf 



u± 3tt 2 (1 + 5) 

(note that the first term does not depend on 5). Then, Eq. (fTo*)) yields the effective potential 
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For 5 = 0, the appearance of combination Wj_(£) + U(£) in potential l|21|l might be expected from the form of the 
last term in Eq. JSJ . Equation (|21|l explicitly shows how the axially nonuniform transverse trapping induces the axial 
potential (which depends on the soliton's "spin"). This result may be realized as transformation of a part of the 
energy of the transverse motion into the axial potential when the confinement frequency gradually varies along x. 

It is well known that, in the notation adopted here, an effective mass of the soliton, if it is treated as a quasiparticle, is 
M e fi = 2N 0| (this relation does not depend on 5). Thus, the equation of motion for the soliton in effective potential 
(|2T|l. M cS (d 2 £/dt 2 ) = -dU cS /dO takes the form 



dt 2 
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This equation clearly shows that, in the absence of the explicit potential, dU/d^ = 0, stable and unstable equilibria 
for the soliton are positions where the local trapping frequency attains its minimum and maximum, respectively. As 
said above, the trapping frequency may also be time-dependent (in addition to being x-dependent), lo±_ = £jj_(£,i), 
then the right-hand side of Eq. (1221 will explicitly depend on t. 

In the most general case, without assuming condition l|19|) to hold, Eq. (|16fl can be solved for a 2 numerically, and 
the effective potential can also be found in a numerical form. To illustrate this possibility, in Fig. 1 we display the 
dependence of a 2 vs. N, as obtained from the numerical solution of Eq. IjKiH at several fixed values of u>±, with 5 = 0. 
As seen from the figure, the transverse width gradually decreases with the increase of N, as predicted by Eq. 120fl . 
The latter effect is a natural manifestation of the self-squeezing of the condensate with attraction between atoms. 

Results obtained by the substitution of the numerically found a 2 in Hamiltonian l|15|l (with 5 = 0), i.e., the effective 
potential, are displayed in Fig. 2. It is seen from this figure that both the linear dependence of the potential on N 
for small N and deviations from the linear dependence at larger N follow the analytical prediction l|21|l . 



IV. CONCLUSIONS 



In this work, we have demonstrated that a tight transverse trap, whose strength gradually varies in the longitudinal 
(axial) direction, induces an effective longitudinal potential for the soliton in the self-attractive BEC. An analytical 
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FIG. 1: The squared transverse size of the condensate, a 2 , as a function of norm N, at several fixed values of the transverse 
trapping frequency, lu±. The curves are obtained from numerical solution of Eq. I|16|l with S = and g 2 = 1, and terminate at 
N = N CI (u>±), as per Eq. JT7J- 




FIG. 2: Dependences of the effective potential 115H with S — 0, U = and = 1 on the soliton's norm N at several fixed 
values of u>±, as obtained by the substitution of numerical data for a 2 , which were used to generate Fig. 1. The plots terminate 
at TV = N CI ( the same way as in Fig. 1. 
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approximation for this potential, based on the variational method, was developed, under the assumption that the 
longitudinal size of the soliton is essentially smaller than the scale of the trap's nonuniformity. A fully explicit result 
was obtained for the case of weak nonlinearity. The potential depends on the soliton's intrinsic vorticity (if any), and 
its norm (the number of atoms trapped in the soliton). The analytical results reported in this Letter can be used 
to analyze inevitable effects of the nonuniformity of the transverse trapping in experimental situations, as well as to 
devise axially nonuniform traps with the purpose to control the longitudinal dynamics of the solitons. Systematic 
numerical simulations of the full 3D model, aimed at direct verification of the analytical predictions, demand special 
effort and will be reported elsewhere. 

It is relevant to mention that a similar problem can also be considered in the model with repulsion (g > in Eq. 
JTJ) and a periodic term if the longitudinal potential generated by an optical lattice, Uo-l(x) — ecos(kx). In that 
case, the system can readily support bright solitons of the gap type 0,0- However, the effective dynamical mass 
of mobile gap solitons in negative [l5j : for this reason, one may expect that, on the contrary to the ordinary solitons 
considered in this Letter, gap solitons will find their stable and unstable equilibria at points where the transverse- 
trapping frequency attains, respectively, its maximum and minimum. These predictions are subject to verification in 
numerical simulations. 
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